• 专利附录
  • 专利说明
  • 权利要求
  • 类似技术
  • 同族专利
  • 引用文献
  • 法律状态
  • 优先权
    • 专利摘要:
    The present disclosure provides a system and method for automatically identifying and classifying seismic terminations within a seismic data volume. A set of surfaces is obtained (step 303) describing the seismic data volume. A plurality of seismic terminations is identified within the set of surfaces (step 307). Based upon seismic attributes or geometric criterion, a termination direction can be determined (step 309) for at least one termination.
    公开号:AU2013219911A1
    申请号:U2013219911
    申请日:2013-01-11
    公开日:2014-07-10
    发明作者:Matthias G. Imhof
    申请人:ExxonMobil Upstream Research Co；
    IPC主号:G01V1-00
    专利说明:
    WO 2013/122704 PCT/US2013/021179 SYSTEM AND METHOD FOR DETECTION AND CLASSIFICATION OF SEISMIC TERMINATIONS CROSS-REFERENCE TO RELATED APPLICATION 5 [0001] This application claims the benefit of U.S. Provisional Patent Application 61/598,233, filed 13 February 2012, entitled SYSTEM AND METHOD FOR DETECTION AND CLASSIFICATION OF SEISMIC TERMINATIONS, the entirety of which is incorporated by reference herein. FIELD OF INVENTION 10 [0002] This invention generally relates to the field of seismic prospecting and reservoir delineation and, more particularly, to a system and method to identify and classify seismic terminations. BACKGROUND [0003] This section is intended to introduce various aspects of the art, which may be 15 associated with exemplary embodiments of the present invention. This discussion is believed to assist in providing a framework to facilitate a better understanding of particular aspects of the present invention. Accordingly, it should be understood that this section should be read in this light, and not necessarily as admissions of prior art. [0004] Seismic stratigraphy analyzes sediments and sedimentary rocks in a geometrical 20 context derived from seismic reflections. A typical first step is the separation of seismic sequence units which is usually done by mapping significant surfaces from where they are indicated by terminations of seismic reflections or seismic reflections with different orientation butting into each other. A second step is the characterization of these seismic sequence units, a process aided by prior determination and classification of the terminations 25 as truncation, toplap, onlap or downlap. Picking and classifying terminations is traditionally a manual task where the interpreter examines the data and marks terminations, for example with a little arrow. [0005] Existing methods for computer-assisted seismic stratigraphy relate to the detection of convergences in seismic images or sections. Representative examples thereof are disclosed 30 by T. Randen et al., "New Seismic Attributes for Automated Stratigraphic Facies Boundary Detection", SEG Technical Program, Expanded Abstracts, pp. 628-631, 1998 and A. Barnes, "Attributes for Automating Seismic Facies Analysis", SEG Technical Program, Expanded Abstracts, pp. 553-556, 2000. - 1- WO 2013/122704 PCT/US2013/021179 [0006] Both of these methods share the need for preliminary determination of a "field of directional vectors" at every pixel of the data. In Barnes' work, quantification of the reflectors convergence is obtained by calculating the divergence of the field of directional vectors. In Randen's work, the directional field is subjected to a search for flowlines whose density 5 expresses the divergence or the convergence of the reflectors. Neither method involves a notion of a surface that separates two stratigraphic units formed under different sedimentological conditions. [0007] In Randen's case, terminations can be defined as locations where flowlines intersect or approximately merge. With an appropriately chosen flowline density, the 10 intersections pinpoint the boundaries by virtually continuous lines, provided sufficient quality and clarity of the seismic data. Randen et al., however, discloses neither further details on requirements for data quality and data clarity nor specifics about flowline generation or intersection detection. In Barnes' case, terminations are best defined in terms of termination density without ability for further classification. 15 [0008] Neither case allows for the automatic detection and classification of individual terminations. Thus, there is a need for improvement in this field. SUMMARY OF THE INVENTION [0009] The present invention provides a system and method for detecting and classifying seismic terminations. 20 [0010] One embodiment of the present disclosure is a computer-implemented method for analyzing a seismic data volume acquired in a seismic survey to determine potential for hydrocarbon accumulations in an associated subsurface region, the method comprising: obtaining surfaces describing the seismic data volume; identifying edge points within the surfaces; identifying a plurality of seismic terminations, wherein each termination is defined 25 by two edge points defining a first terminated surface and a second terminated surface; and using a computer to determine a termination direction for at least one termination. [0011] Another embodiment of the present disclosure is a computer-implemented method for analyzing a seismic data volume acquired in a seismic survey to determine potential for hydrocarbon accumulations in an associated subsurface region, the method comprising: 30 obtaining a set of surfaces describing the seismic data volume; identifying a plurality of seismic terminations within the set of surfaces; and using a computer to determine a termination direction for at least one termination. -2- WO 2013/122704 PCT/US2013/021179 [0012] The foregoing has broadly outlined the features of some embodiments of the present invention in order that the detailed description that follows may be better understood. Additional features and embodiments will also be described herein. BRIEF DESCRIPTION OF THE DRAWINGS 5 [0013] The present invention and its advantages will be better understood by referring to the following detailed description and the attached drawings. [0014] Figures 1A-1D depict examples of the four classic stratigraphic termination types: truncation (Figure 1A), toplap (Figure 1B), onlap (Figure IC), and downlap (Figure ID). [0015] Figure 2 depicts an example of a termination of two seismic surfaces. 10 [0016] Figure 3 is a flow chart showing the basic steps of automatic termination type classification according to one embodiment of the present disclosure. [0017] Figure 4 demonstrates the recognition of an upward termination according to one embodiment of the present disclosure. [0018] Figure 5 demonstrates the recognition of a downward termination according to 15 one embodiment of the present disclosure. [0019] Figure 6 demonstrates the recognition of an undefined termination direction according to one embodiment of the present disclosure. [0020] Figures 7A-7D demonstrate the use of relative steepness of terminated and persistent surfaces to differentiate between truncation (Figure 7A), toplap (Figure 7B), onlap 20 (Figure 7C), and downlap (Figure 7D). [0021] Figure 8 depicts a seismic cross section from a processed seismic data volume. [0022] Figure 9 depicts the seismic cross section of Figure 8 superimposed with truncations which have been identified according to one embodiment of the present disclosure. 25 [0023] Figure 10 depicts the seismic cross section of Figure 8 superimposed with toplaps which have been identified according to one embodiment of the present disclosure. [0024] Figure 11 depicts the seismic cross section of Figure 8 superimposed with onlaps which have been identified according to one embodiment of the present disclosure. [0025] Figure 12 depicts the seismic cross section of Figure 8 superimposed with 30 downlaps which have been identified according to one embodiment of the present disclosure. -3 - WO 2013/122704 PCT/US2013/021179 [0026] Figure 13 depicts the seismic cross section of Figure 8 superimposed with upward terminations which have been identified according to one embodiment of the present disclosure. [0027] Figure 14 depicts the seismic cross section of Figure 8 superimposed with 5 downward terminations which have been identified according to one embodiment of the present disclosure. [0028] Figure 15 depicts a truncation density section which has been determined according to one embodiment of the present disclosure. [0029] It should be noted that the figures are merely examples of several embodiments of 10 the present invention and no limitations on the scope of the present invention are intended thereby. Further, the figures are generally not drawn to scale, but are drafted for purposes of convenience and clarity in illustrating various aspects of certain embodiments of the invention. DESCRIPTION OF THE SELECTED EMBODIMENTS 15 [0030] For the purpose of promoting an understanding of the principles of the invention, reference will now be made to the embodiments illustrated in the drawings and specific language will be used to describe the same. It will nevertheless be understood that no limitation of the scope of the invention is thereby intended. Any alterations and further modifications in the described embodiments, and any further applications of the principles of 20 the invention as described herein are contemplated as would normally occur to one skilled in the art to which the invention relates. One embodiment of the invention is shown in great detail, although it will be apparent to those skilled in the relevant art that some features that are not relevant to the present invention may not be shown for the sake of clarity. [0031] Persons skilled in the technical field will readily recognize that in practical 25 applications of the disclosed methodology, it must be performed on a computer, typically a suitably programmed digital computer. Further, some portions of the detailed descriptions which follow are presented in terms of procedures, steps, logic blocks, processing and other symbolic representations of operations on data bits within a computer memory. These descriptions and representations are the means used by those skilled in the data processing 30 arts to most effectively convey the substance of their work to others skilled in the art. In the present application, a procedure, step, logic block, process, or the like, is conceived to be a self-consistent sequence of steps or instructions leading to a desired result. The steps -4- WO 2013/122704 PCT/US2013/021179 sometimes require physical manipulations of physical quantities. Usually, although not necessarily, these quantities take the form of electrical or magnetic signals capable of being stored, transferred, combined, compared, and otherwise manipulated in a computer system. [0032] It should be borne in mind, however, that all of these and similar terms are to be 5 associated with the appropriate physical quantities and are merely convenient labels applied to these quantities. Unless specifically stated otherwise as apparent from the following discussions, it is appreciated that throughout the present application, discussions utilizing the terms such as "obtaining", "generating", "identifying", "determining", and "displaying" or the like, refer to the action and processes of a computer system, or similar electronic 10 computing device, that manipulates and transforms data represented as physical (electronic) quantities within the computer system's registers and memories into other data similarly represented as physical quantities within the computer system memories or registers or other such information storage, transmission or display devices. [0033] Seismic stratigraphy is a branch of stratigraphy in which sediments and 15 sedimentary rocks are interpreted in a geometrical context from seismic reflectors. Seismic sequence analysis separates seismic data into seismic-sequence units or seismic-facies units. This separation is often done by mapping unconformities where they are indicated by angularity or seismic reflections with different orientation terminating each other. The unconformities are then extended to where they cannot be identified in this manner by 20 following along reflections, advantage being taken of the fact that the unconformity reflection is often relatively strong. The interpretation of the seismic-sequence units is aided by analysis of the terminations, or locations where reflections end. [0034] In order to aid in the understanding of the present disclosure, a brief discussion of terminations is provided below. As understood by those of ordinary skill in the art, a 25 termination type is based on the manner in which reflections come to an end. Lapout is the lateral termination of a reflection or associated stratal unit when it becomes too thin to be resolved. Truncation, on the other hand, is the lateral termination of a reflection or associated stratal unit when other strata cut it off. [0035] Four classic termination types are depicted in Figures 1A-ID. Figure 1A depicts a 30 truncation, Figure lB depicts a toplap, Figure IC depicts an onlap, and Figure ID depicts a downlap. As demonstrated by the figures, truncations and toplaps are terminations against an upper boundary, or unit, 101, while onlap and downlap are terminations against a lower boundary 103. -5 - WO 2013/122704 PCT/US2013/021179 [0036] As demonstrated in Figure 1A, truncation is the lateral termination of a reflection by erosion. It occurs at the upper boundary 101 of a depositional sequence and may extend over a wide regional area or be confined to a local scour. As used throughout the present disclosure, a truncation is defined as termination of a relatively flat reflection by a relatively 5 steep reflection. [0037] Structural truncation is the lateral termination of a unit, or surface, by structural disruption. Such truncation is most easily recognized where it cuts across a stack of concordant units or surfaces. Often, the disruption is caused by faulting. At times, the distinction between erosional and structural truncation can be difficult. For the purpose of 10 understanding the principles of the disclosed system and method, erosional and structural truncations are not distinguished and may both be referred to as a truncation. [0038] As demonstrated in Figure IB, a toplap is a lapout at the upper boundary 101 of a depositional unit. Initially inclined strata, such as clinoforms or bar forms, may show this relation. Often, lateral terminations may taper, updip and asymptotically approach the upper 15 boundary. Due to limited resolution on seismic sections, reflections may appear to terminate abruptly at a high angle against the upper, relatively horizontal surface. As used throughout the present disclosure, a toplap is defined as the updip termination of an inclined reflection against a relatively flat reflection. [0039] Subsequent structural deformation can change the reflection orientations and may 20 even convert a truncation to an apparent toplap or may convert a toplap to an apparent truncation. As used throughout the present disclosure, the term "upward termination" is used to describe a termination against an upper unit or surface 101. Thus, an upward termination denotes either an (apparent) truncation or an (apparent) toplap. Throughout this disclosure, the word "apparent" is typically suppressed because all termination terminations and 25 classifications are apparent since they are based primarily on geometry, dip and/or other seismic attributes. [0040] For the case of an onlap and as depicted in Figure 1 C, either an initially horizontal unit laps out against an initially inclined surface, or an initially inclined unit laps out updip against a surface of greater initial inclination. For the case of downlap and as depicted in 30 Figure ID, an initially inclined unit terminates downdip against an initially horizontal surface or surface of lesser initial inclination. Onlap or downlap can usually be readily identified based on their differences in orientation. Subsequent structural deformation can change their orientations and may even convert an onlap to an apparent downlap or may convert a -6- WO 2013/122704 PCT/US2013/021179 downlap to an apparent onlap. As used throughout the present disclosure, the term "downward termination" is used to describe a termination against a lower unit or surface. Thus, a downward termination denotes either a downlap or an onlap. [0041] As used throughout the present disclosure, a seismic termination is defined as the 5 location where a seismic event such as a peak (local maximum) or trough (local minimum) ends. In a seismic cross section or seismic volume, however, an event rarely just ends. Rather, the event terminates or disappears by merging with another event. However, a seismic peak being a local maximum on a seismic trace cannot disappear unless a vertically adjacent trough disappears, too. Otherwise, there would be two vertically adjacent local 10 minima left without an in-between local maximum, which is impossible. Similarly, a seismic trough being a local minimum on a seismic trace cannot disappear unless a vertically adjacent peak or local maximum disappears, too. Otherwise, there would be two vertically adjacent local maxima left without an in-between local minimum, which is again impossible. Thus, events defined by correlated peaks and troughs terminate in pairs. 15 [0042] Figure 2 illustrates this behavior with a synthetic seismic wedge model. Persistent surfaces defined by correlated trough events are designated with light dashed lines 201 and persistent surfaces defined by correlated peak events are designated with light solid lines 203. The adjacent surfaces drawn as heavy dashed lines 205 and solid lines 207 (correlated trough or peak events, respectively) terminate at the same lateral location. The edge points of 20 surfaces 205 and 207 are located in adjacent peak and trough events of the same trace, and thus are deemed to be valid terminations. If the surfaces 205 and 207 had been mistracked by splitting them into smaller patches, then the edge points are likely to be non-adjacent and/or in different traces which allows classifying these edge points as being invalid. As will be discussed in more detail below, one aspect of the present disclosure is the determination of 25 whether 205 and 207 are truncated by the light solid line 209 or by the light dashed line 211, i.e., determination whether the termination direction is upwards or downwards. It is easy to see that 205 and 207 merge into 209 while 211 does not show any indication on the presence and termination of 205 and 207. Thus, the termination direction is determined to be upwards, surface 209 is termed the persistent surface, and surfaces 205 and 207 are termed the 30 terminated surfaces. [0043] The flow chart of Figure 3 will be referred to in describing one embodiment of the present disclosure. While each step will be described in further detail below, the depicted process (301) first obtains a set of surfaces (303) which describe a seismic data volume -7- WO 2013/122704 PCT/US2013/021179 acquired in a seismic survey. Edge points are then identified (305) within the set of surfaces. Valid terminations are identified (307) and the valid terminations are then evaluated to determine a termination direction (309). The process may continue by further evaluating the terminations to determine and assign a termination type (311) to each valid termination. A 5 user may then select at least one particular termination direction or type (313). The process may continue by determining a termination density (315) for the selected termination direction or type. The selected termination direction, type, and/or density may then be displayed to the user (317). In some embodiments of the present disclosure, the order of steps 303-317 may be changed, at least one step may be repeated, and/or at least one step 10 may not be performed. For example, one embodiment of the present disclosure performs only steps 303-307. [0044] As noted above, one embodiment of the disclosed method begins by obtaining a set of surfaces (303). U.S. Patent Application No. 12/920,013, titled "Seismic Horizon Skeletonization" to Imhof et al., which is incorporated by reference in its entirety, discloses 15 an algorithm that extracts essentially every event from a seismic dataset and assigns these events to surfaces. In one embodiment, such skeletonization techniques are used to derive seismic surfaces from a seismic data volume. Different embodiments may extract or obtain seismic surfaces using other techniques, such as, but not limited to, tracking peaks and/or troughs or either kind of zero crossings (+/- or -/+). Other methods for generating or 20 obtaining a large set of surfaces from a seismic data volume include those taught by U.S. Pat. No. 7,248,539, titled "Extrema Classification" to Borgos; U.S. Pat. No. 5,570,106, titled "Method and Apparatus for Creating Horizons from 3-D Seismic Data" to Viswanathan; and U.S. Patent Application Publication No. 2008/0285384, titled "System and Method for Displaying Seismic Horizons with Attributes" to James. The disclosed system and method is 25 not limited to a particular way of tracking, generating or obtaining surfaces. In one embodiment, all peak and trough surfaces are extracted. [0045] At step 305, edge points are then identified within the set of surfaces. Depending on the desired connectivity of the surface, an edge point may be defined as a point that is not surrounded by four or eight points that belong to the same surface. With this edge point 30 definition, every tiny hole inside a given surface creates additional edge points. In one embodiment, holes within surfaces are left untouched as many holes attributed to mistracked events will be deemed invalid edge points because they violate the requirement of two adjacent edge points within one trace. Small perturbations in amplitudes and waveforms can -8- WO 2013/122704 PCT/US2013/021179 also lead to holes within surfaces. In some embodiments, interior holes are filled and thus removed, for example, by interpolation. In some embodiments, only a subset of the interior holes are filled, such as, but not limited to, holes that match a user-specified criteria. Such user-specified criteria may include holes having an insufficient size and/or having a particular 5 shape. [0046] Valid terminations are identified at step 307. In one embodiment, valid terminations are determined by filtering edge points determined to be unassociated with stratigraphic terminations. One potential filtering operation is to discard edge points interior to a surface. Another potential filtering operation is to remove edge points that do not occur 10 strictly in vertically-stacked pairs. Figure 2 demonstrates that acceptable terminations occur in pairs because both a peak and a trough surface terminate at the same location. In the presence of noise or coherent imaging artifacts, the skeletonization algorithm disclosed in U.S. Application No. 12/920,013 can break surfaces into smaller pieces, thereby generating artificial terminations along their edges. Such artificial terminations can be recognized 15 because there rarely is one vertically adjacent companion termination as required for a valid termination. Another filtering operation relates to seismic noise that can introduce vertical trains of events, and thus trains of potential terminations. Nevertheless, such noise-train terminations can be distinguished from valid terminations because they appear in vertical trains and not just in simple pairs as required for valid stratigraphic terminations. In some 20 embodiments, a majority of the edge points and potential terminations identified within the obtained set of surfaces can be discarded as a result of step 307 because they do not appear in strictly vertically adjacent pairs. Regardless of the filtering operation utilized, the remaining edge points are, therefore, associated with valid stratigraphic terminations. [0047] Process 301 proceeds by determining a termination direction (step 309) for each 25 valid termination. In other words, the next step after the identification of valid edge points and terminations (or after the rejection of invalid edge points and terminations) is the recognition of whether a terminating event merges with the upper or lower event. In some embodiments, the termination direction is determined for some, but not all, valid terminations. The determination of termination direction may be based on a seismic attribute 30 and/or geometric criterion. In some embodiments, seismic amplitude is used to determine the termination direction by directing the termination toward the stronger of the two adjacent persistent surfaces. In some embodiments, the termination is directed toward the adjacent persistent surface that exhibits a broader waveform or toward the one that exhibits a lower -9- WO 2013/122704 PCT/US2013/021179 dominant frequency. In some embodiments, a waveform attribute is used to direct the termination. One example of such an embodiment exploits an anomaly of the waveform to assign the termination direction, namely the presence of a doubleted waveform. [0048] U.S. Patent Application 12/920,285, titled "Method For Geophysical And 5 Stratigraphic Interpretation Using Waveform Anomalies", to Imhof, which is incorporated by reference in its entirety, discloses an interpretation method that exploits an anomalous waveform to determine locations where a seismic wavelet separates to create a set of two additional events. At times, a doubleted waveform appears when one seismic wavelet is about to separate into two events. Such a doubleted waveform constitutes a waveform anomaly. In 10 one embodiment, the location and/or geometry of the wavelet separation is used to determine termination direction. [0049] Figure 4 shows a vanishing waveform-maximum event 401 terminating in an ordinary maximum. Vanishing waveform-minimum event 403, however, ends in a doublet. As used throughout the present disclosure, a doublet is either a local waveform minimum 15 with positive value or a local waveform maximum with a negative value. In contrast, a nondoublet is used herein to describe an ordinary maximum or minimum event. [0050] With respect to Figure 4, event 401 is considered a nondoublet. In this case, event 403 is a local waveform minimum with a positive value, and thus a doublet. As illustrated, events 401 and 403 are likely to be terminated by the above waveform-maximum event 405 20 and not the beneath waveform-minimum event 407. In other words, events 401, 403 and 405 all feed into the same wiggle at the termination, and events 401 and 403 are likely to terminate into 405. On the other hand, events 401, 403, and 407 feed into two different wiggles of opposing polarity; therefore, event 401 and 403 are unlikely to terminate into 407. Thus, the upward termination 409 is placed where events 401 and 403 vanish. Note that 25 changing the waveform polarities also leads to an upward termination. [0051] Figure 5 shows the case of a vanishing waveform-maximum event 501 terminating in an ordinary maximum. Vanishing waveform minimum 503, however, ends in a doublet of a local waveform-minimum with positive value. Due to the arrangement of the terminated surfaces, events 501 and 503 are likely to be terminated by the beneath waveform 30 maximum event 505 and not the above waveform-minimum event 507. Thus, downward termination 509 is placed where events 501 and 503 vanish. Again, note that changing the waveform polarities leads to a downward termination. - 10 - WO 2013/122704 PCT/US2013/021179 [0052] Therefore, in one embodiment of the present disclosure, the presence of the doublet waveform allows determination of the termination direction. For two vertically adjacent edge points where one of the edge points lies on a waveform doublet and the other edge point lies on a waveform nondoublet, the termination direction is from the nondoublet to 5 the doublet. Thus, if the nondoublet is below the doublet, then the termination direction is upward. However, if the nondoublet is above the doublet, then the termination direction is downward. [0053] Figure 6 shows a case where the waveform-minimum event 601 ends with a regular minimum and the waveform-maximum event 603 ends with a regular maximum. 10 Thus, the termination direction cannot be determined by exploitation of the doublet waveform. In one embodiment, the termination is deemed to be an undefined termination with undefined direction. In another embodiment, event 601 may be deemed an upward termination, while event 603 is deemed a downward termination. As a result, the termination defined by events 601 and 603 is defined as both upward and downward. Lacking a 15 doubleted waveform to determine the termination direction, some embodiments of the present disclosure may employ other properties of the obtained surfaces or other seismic attributes to determine termination direction. In some embodiments, the user may elect to use another criterion even if doubleted waveforms are available. In one embodiment, events 601 and 603 may be deemed to be truncated by the stronger of events 605 or 607. An event may be found 20 to be "stronger" based on its seismic amplitude within the seismic trace containing the terminated surface edge points. In another embodiment, an event may be found to be "stronger" based on its seismic amplitudes average for events belonging to the same surface but within a user-specified distance from the termination. A particular case of this embodiment is to define event strength by forming the average of all events belonging to the 25 same surface. Another example is the use of seismic loop duration where the termination direction for the 601 and 603 events are defined toward the surface 605 or 607 with the broader wavelet, i.e., the wavelet with the longer duration or lower dominant frequency. Indeed, any attribute, such as, but not limited to, amplitude, shape, or geometry, may be used to define the direction of termination. 30 [0054] Returning again to process 301, step 311 determines the termination type for at least one termination. More specifically, the upward terminations may be classified as truncations or toplaps, while the downward terminations may be classified as onlaps or downlaps. In one embodiment of the present disclosure, the termination type determination - 11 - WO 2013/122704 PCT/US2013/021179 is made based on a comparison of the relative surface orientations, i.e., the two surfaces that are terminated and the persistent surface that terminates the others. In such an embodiment, the first step is to estimate the surface normal direction of at least one of the terminated surfaces as well as the persistent surface near the termination. The normal directions may be 5 determined by fitting the nearby points with a plane. The directions of these surfaces are then compared. Based on this comparison as well as the termination direction, the termination type may then be determined and assigned. [0055] Figures 7A-7D demonstrate the use of relative steepness of terminated and persistent surfaces to differentiate termination type. Figure 7A shows surface 701 cutting off 10 surfaces 702 and 703 from above. Because the dip 704 of surface 701 is greater than the dip 705 of surfaces 702 and 703, the termination 700 is classified to be an truncation. Therefore, for a truncation, the persistent surface is steeper than the terminated surfaces. [0056] Figure 7B shows surface 711 cutting off surfaces 712 and 713 from above. Because the dip 714 of surface 711 is lesser than the dip 715 of surfaces 712 and 713, the 15 termination 710 is classified to be a toplap. Therefore, for toplaps, the terminated surfaces are steeper than the persistent one. In Figure 7C, surface 721 cuts off surfaces 722 and 723 from below. Because the dip 724 of surface 721 is greater than the dip 725 of surfaces 722 and 723, the termination 720 is classified to be an onlap. Therefore, for onlaps, the persistent surface is steeper than the terminated ones. Lastly and as shown in Figure 7D, surface 731 20 cuts off surfaces 732 and 733 from below. Because the dip 734 of surface 731 is lesser than the dip 735 of surfaces 732 and 733, the termination 730 is classified to be a downlap. Therefore, for downlaps, the terminated surfaces are steeper than the persistent one. [0057] Proceeding through process 301 depicted in Figure 3, a user then selects at least one termination direction or type (step 313) of interest, such as upward terminations or 25 toplaps and onlaps. In some embodiments, the user also specifies a window size. In one embodiment, the user may input his/her selection into a computer or other electronic device using user interfaces known to or understood by those of ordinary skill in the art. [0058] At step 315, a termination density is determined for the selected termination direction or type. In one embodiment, an algorithm counts at specified points how many 30 terminations of the specified direction or type are contained in a window of the specified size anchored at that point. In some embodiments, essentially every point of the seismic data volume anchors such a window. In some embodiments, the count is unweighted thereby providing the actual number of terminations of the specified type within the window. In - 12 - WO 2013/122704 PCT/US2013/021179 another embodiment, the contribution of each termination is weighted, such as, but not limited to, with respect to its distance from the window center or anchor point. Either embodiment may be implemented by convolution of a binary indicator volume for the specified termination direction or type with an appropriate filter. 5 [0059] Step 317 displays the results for review or analysis. In one embodiment, terminations of the selected termination type are superimposed over a seismic section extracted from the seismic data volume. One suitable method of visualization is optical smashing of the individual terminations. In an optical smash, terminations outside a specified region are rendered transparent. Terminations within that region are rendered and typically 10 color coded by termination type. Optical smashing demonstrates that terminations tend to cluster by termination type. In one embodiment, the optically smashed terminations are co rendered with seismic amplitudes or a seismic section extraction from the seismic data volume, which facilitates high-resolution seismo-stratigraphic interpretation of essentially every seismic surface. 15 [0060] As described herein, terminations are located and defined on the terminated surfaces. In one embodiment of the present disclosure, the terminations are relocated from the terminated surfaces to the persistent surfaces which cut the terminated surfaces. Relocating the terminations in this manner allows classifying surfaces as unconformities or flooding surfaces. In one embodiment, a surface that predominantly contains relocated truncations is 20 classified as an unconformity. A surface that predominantly contains relocated downlaps is likely to be a flooding surface. [0061] Figure 8 shows an example of a seismic section extracted from a larger three dimensional seismic data volume. In this example, the seismic surfaces were computed by application of a skeletonization algorithm to the entire seismic data volume. The resulting 25 surfaces are evaluated pursuant to an embodiment of the disclosed method to determine and classify seismic terminations. Figure 9 shows the same cross section but with truncations from the nine nearest slices projected, or superimposed, onto the cross section. In Figure 9, the depicted truncations represent an optical smash over nine seismic data slices. Figure 10 shows an optical smash of the identified toplaps. Figure 11 depicts an optical smash of the 30 identified onlaps, whereas Figure 12 depicts an optical smash of the identified downlaps. Figures 13 and 14 present upward and downward terminations, respectively. The classification into upward and downward terminations is more robust than the classification into truncations, toplaps, onlaps and downlaps because they do not require estimation and - 13 - WO 2013/122704 PCT/US2013/021179 comparison of local dips. Figure 15 presents the truncation density obtained by counting the number of truncations within probes with a size of 21 by 21 by 7 voxels. [0062] In one embodiment of the present disclosure and as described above, a set of surfaces is obtained and edge points are then identified within that set of surfaces. In other 5 embodiments, anomalies in the obtained seismic waveforms are used to seed an automatic event tracker. The automatic event tracker is used to create a set of surfaces which are anchored at the location of the anomalous waveform. Based on the principles described herein, these surfaces may then be analyzed for termination direction and/or termination type. As presented herein, waveform doublets may be used to determine termination direction. In 10 some embodiments, terminations not located within doubleted waveforms are left with an undefined termination direction. In an alternative embodiment, the seismic data is first scanned for doubleted waveforms. In such an embodiment, instead of tracking essentially every surface in the data, only those surfaces that touch a doubleted waveform are tracked. Said differently, the doubleted waveforms are utilized to seed the tracker. 15 [0063] Embodiments of the present invention also relate to an apparatus for performing the operations herein. This apparatus may be specially constructed for the required purposes, or it may comprise a general-purpose computer selectively activated or reconfigured by a computer program stored in the computer. Such a computer program may be stored in a computer readable medium. A computer-readable medium includes any mechanism for 20 storing or transmitting information in a form readable by a machine (e.g., a computer). For example, but not limited to, a computer-readable (e.g., machine-readable) medium includes a machine (e.g., a computer) readable storage medium (e.g., read only memory ("ROM"), random access memory ("RAM"), magnetic disk storage media, optical storage media, flash memory devices, etc.), and a machine (e.g., computer) readable transmission medium 25 (electrical, optical, acoustical or other form of propagated signals (e.g., carrier waves, infrared signals, digital signals, etc.)). [0064] Furthermore, as will be apparent to one of ordinary skill in the relevant art, the modules, features, attributes, methodologies, and other aspects of the invention can be implemented as software, hardware, firmware or any combination of the three. Of course, 30 wherever a component of the present invention is implemented as software, the component can be implemented as a standalone program, as part of a larger program, as a plurality of separate programs, as a statically or dynamically linked library, as a kernel loadable module, as a device driver, and/or in every and any other way known now or in the future to those of - 14 - WO 2013/122704 PCT/US2013/021179 skill in the art of computer programming. Additionally, the present invention is in no way limited to implementation in any specific operating system or environment. [0065] It should be understood that the preceding is merely a detailed description of specific embodiments of this invention and that numerous changes, modifications, and 5 alternatives to the disclosed embodiments can be made in accordance with the disclosure here without departing from the scope of the invention. The preceding description, therefore, is not meant to limit the scope of the invention. Rather, the scope of the invention is to be determined only by the appended claims and their equivalents. It is also contemplated that structures and features embodied in the present examples can be altered, rearranged, 10 substituted, deleted, duplicated, combined, or added to each other. The articles "the", "a" and "an" are not necessarily limited to mean only one, but rather are inclusive and open ended so as to include, optionally, multiple such elements. 15 - 15 -
    权利要求:
    Claims (31)
    [1] 1. A computer-implemented method for analyzing a seismic data volume acquired in a seismic survey to determine potential for hydrocarbon accumulations in an associated subsurface region, the method comprising: 5 obtaining surfaces describing the seismic data volume; identifying edge points within the surfaces; identifying a plurality of seismic terminations, wherein each termination is defined by two edge points defining a first terminated surface and a second terminated surface; and using a computer to determine a termination direction for at least one termination. 10
    [2] 2. The method of claim 1, wherein the obtained surfaces are derived from the seismic data volume using skeletonization.
    [3] 3. The method of claim 1, wherein the obtained surfaces are derived from tracking peaks 15 and troughs within the seismic data volume.
    [4] 4. The method of claim 1 further comprising scanning the seismic data volume for doubleted waveforms, wherein the obtained surfaces seeded at the doubleted waveforms. 20
    [5] 5. The method of claim 1, wherein the first terminated surface ends at a doublet and the second terminated surface ends at a nondoublet, the termination direction is defined in the direction from the second surface to the first surface.
    [6] 6. The method of claim 1, wherein the termination direction is selected from a group 25 consisting of upward, downward, and undefined.
    [7] 7. The method of claim 1 further comprising determining a termination type to the at least one termination. 30
    [8] 8. The method of claim 7, wherein the first and second terminated surfaces each end on or adjacent to a persistent surface
    [9] 9. The method of claim 8, wherein determining the termination type is based upon the geometric orientations of the first terminated surface and the persistent surface. - 16 - WO 2013/122704 PCT/US2013/021179
    [10] 10. The method of claim 8, wherein the termination type is selected from a group consisting of truncation, toplap, onlap, and downlap. 5
    [11] 11. The method of claim 9, wherein determining the termination type comprises: estimating a first angle defined by a surface normal direction of the first terminated surface; and estimating a second angle defined by a surface normal direction of the persistent surface, 10 wherein the termination type is based upon the termination direction, the first angle, and the second angle.
    [12] 12. The method of claim 10 further comprising: relocating the terminations from the end of the terminate surfaces to the associated 15 persistent surface; and classifying the persistent surface as an unconformity or a flooding surface based on the relocated termination types defined along the persistent surface.
    [13] 13. The method of claim 12, wherein the persistent surface predominately containing 20 relocated truncations is classified as an unconformity.
    [14] 14. The method of claim 12, wherein the persistent surface predominately containing relocated downlaps is classified as a flooding surface. 25
    [15] 15. The method of claim 1 further comprising selecting a termination type.
    [16] 16. The method of claim 15 further comprising displaying at least one termination of the selected termination type superimposed over a seismic section extracted from the seismic data volume. 30
    [17] 17. The method of claim 15 further comprising determining a termination density of the selected termination type.
    [18] 18. A method for producing hydrocarbons comprising: - 17 - WO 2013/122704 PCT/US2013/021179 conducting a seismic survey of a subsurface region; obtaining an interpretation of the subsurface using the method of claim 1; drilling a well into the subsurface based at least in part on the interpretation; and producing hydrocarbons from the well. 5
    [19] 19. A computer-implemented method for analyzing a seismic data volume acquired in a seismic survey to determine potential for hydrocarbon accumulations in an associated subsurface region, the method comprising: obtaining a set of surfaces describing the seismic data volume; 10 identifying a plurality of seismic terminations within the set of surfaces; and using a computer to determine a termination direction for at least one termination.
    [20] 20. The method of claim 19, further comprising identifying edge points within the set of surfaces, wherein each termination is defined by two edge points defining a first terminated 15 surface and a second terminated surface.
    [21] 21. The method of claim 20, wherein the first terminated surface ends at a doublet and the second terminated surface ends at a nondoublet, the termination direction is defined in the direction from the second surface to the first surface. 20
    [22] 22. The method of claim 20, wherein the termination direction is selected from a group consisting of upward, downward, and undefined.
    [23] 23. The method of claim 19 further comprising determining a termination type to the at 25 least one termination.
    [24] 24. The method of claim 23, wherein the first and second terminated surfaces each end on or adjacent to a persistent surface 30
    [25] 25. The method of claim 23, wherein determining the termination type is based upon the geometric orientations of the first terminated surface and the persistent surface.
    [26] 26. The method of claim 23, wherein the termination type is selected from a group consisting of truncation, toplap, onlap, and downlap. 35 - 18 - WO 2013/122704 PCT/US2013/021179
    [27] 27. The method of claim 23, wherein determining the termination type comprises: estimating a first angle defined by a surface normal direction of the first terminated surface; and estimating a second angle defined by a surface normal direction of the persistent 5 surface, wherein the termination type is based upon the termination direction, the first angle, and the second angle.
    [28] 28. The method of claim 1, wherein each termination is defined by two vertically adjacent 10 edge points defining a first terminated surface and a second terminated surface.
    [29] 29. The method of claim 19, further comprising identifying edge points within the set of surfaces, wherein each termination is defined by two vertically adjacent edge points defining a first terminated surface and a second terminated surface. 15
    [30] 30. The method of claim 8, wherein determining the termination type is based upon geometric orientations of the persistent surface and at least one of the first terminated surface or the second terminated surface. 20
    [31] 31. The method of claim 30, wherein determining the termination type comprises: estimating a first angle defined by a surface normal direction of the first terminated surface, a surface normal direction of the second terminated surface, or the surface normal direction of the first and second terminated surfaces; and 25 estimating a second angle defined by a surface normal direction of the persistent surface, wherein the termination type is based upon the termination direction, the first angle, and the second angle. 30 - 19 -
    类似技术:
    公开号 | 公开日 | 专利标题
    AU2013219911B2|2014-11-06|System and method for detection and classification of seismic terminations
    US8903659B2|2014-12-02|Method for analysis seismic data with attribute functions
    US9366772B2|2016-06-14|Method for creating a hierarchically layered earth model
    US20100214870A1|2010-08-26|Method and apparatus for dynamic extraction of extrema-based geometric primitives in 3d voxel volumes
    US7257488B2|2007-08-14|Method of sedimentologic interpretation by estimation of various chronological scenarios of sedimentary layers deposition
    US20130138350A1|2013-05-30|Enhancedtermination identification function based on dip field generated from surface data
    Basir et al.2013|Multi-attribute ant-tracking and neural network for fault detection: A case study of an Iranian oilfield
    CN108226999A|2018-06-29|The processing method of the small scale fracture hole body information of carbonate rock
    US20090192717A1|2009-07-30|Methods and apparatuses for dynamic filtering of geometric primitives in 3d space
    CN104011566A|2014-08-27|Wavelet-transform based system and method for analyzing characteristics of a geological formation
    Lacaze et al.2011|Seismic stratigraphic interpretation from a geological model-A north sea case study
    US20100228486A1|2010-09-09|Method and system for seismic data processing
    CN104781697A|2015-07-15|System and method for velocity anomaly analysis
    Pelissier et al.2017|Interpretation value of diffractions and sub-specular reflections–applications on the Zhao Dong field
    Sicking et al.2014|Case study comparing micro-earthquakes, fracture volumes, and seismic attributes
    Refayee et al.2016|Fault and fracture detection in unconventional reservoirs: A Utica shale study
    Rush et al.2017|Geostatistical facies modeling trends for oolitic tidal sand shoals
    Maroquin et al.2009|A visual data-mining methodology to conduct seismic facies analysis: Part 2—Testing and comparison with other unsupervised clustering methods
    Hammon III2009|Computer techniques to aid the interpretation of salt bodies and stratigraphy in three-dimensional seismic volumes
    Wetterhahn2014|Fracture geometry and orientation correlations in a fold-and-thrust structure in the Northern Apennines, Italy
    Johnson2019|Utilizing Geophysical Attributes to Investigate the Architecture of a Pinnacle Reef Complex, Michigan Basin, USA
    Sodagar et al.2010|2D seismic modeling of the Redwater Leduc Reef, Alberta
    Liu et al.2021|Convolutional neural network for guided-wave energy identification in microseismic DAS data
    Boschetti et al.2000|Improved edge detection in potential field maps and graphical estimation of depth-to-the-top
    Pate2012|Textural analysis of FMI logs to determine lithofacies in the early paddock member
    同族专利:
    公开号 | 公开日
    US20140345857A1|2014-11-27|
    CA2861830A1|2013-08-22|
    WO2013122704A1|2013-08-22|
    US9523782B2|2016-12-20|
    CA2861830C|2015-05-05|
    EP2815255B1|2017-03-01|
    EP2815255A4|2015-10-14|
    EP2815255A1|2014-12-24|
    AU2013219911B2|2014-11-06|
    引用文献:
    公开号 | 申请日 | 公开日 | 申请人 | 专利标题
    US4916615A|1986-07-14|1990-04-10|Conoco Inc.|Method for stratigraphic correlation and reflection character analysis of setsmic signals|
    US5047991A|1989-04-28|1991-09-10|Schlumberger Technology Corporation|Lithology identification using sonic data|
    US4992995A|1989-10-24|1991-02-12|Amoco Corporation|Methods for attenuating noise in seismic data|
    US5274714A|1990-06-04|1993-12-28|Neuristics, Inc.|Method and apparatus for determining and organizing feature vectors for neural network recognition|
    US5265192A|1990-09-20|1993-11-23|Atlantic Richfield Company|Method for the automated editing of seismic traces using an adaptive network|
    US5570106A|1993-07-28|1996-10-29|Landmark Graphics Corporation|Method and apparatus for creating horizons from 3-D seismic data|
    US5444619A|1993-09-27|1995-08-22|Schlumberger Technology Corporation|System and method of predicting reservoir properties|
    US5416750A|1994-03-25|1995-05-16|Western Atlas International, Inc.|Bayesian sequential indicator simulation of lithology from seismic data|
    GB2293010B|1994-07-07|1998-12-09|Geco As|Method of processing seismic data|
    US5586082A|1995-03-02|1996-12-17|The Trustees Of Columbia University In The City Of New York|Method for identifying subsurface fluid migration and drainage pathways in and among oil and gas reservoirs using 3-D and 4-D seismic imaging|
    US5539704A|1995-06-23|1996-07-23|Western Atlas International, Inc.|Bayesian sequential Gaussian simulation of lithology with non-linear data|
    FR2738920B1|1995-09-19|1997-11-14|Elf Aquitaine|METHOD FOR AUTOMATIC SEISMIC FACIAL RECOGNITION|
    US6012017A|1996-09-25|2000-01-04|Geoquest, A Division Of Schlumberger Technology Corporation|Interpreting formation tops|
    US6052650A|1997-02-27|2000-04-18|Schlumberger Technology Corporation|Enforcing consistency in geoscience models|
    US6466923B1|1997-05-12|2002-10-15|Chroma Graphics, Inc.|Method and apparatus for biomathematical pattern recognition|
    GB9904101D0|1998-06-09|1999-04-14|Geco As|Subsurface structure identification method|
    GB9819910D0|1998-09-11|1998-11-04|Norske Stats Oljeselskap|Method of seismic signal processing|
    US6574565B1|1998-09-15|2003-06-03|Ronald R. Bush|System and method for enhanced hydrocarbon recovery|
    US6236942B1|1998-09-15|2001-05-22|Scientific Prediction Incorporated|System and method for delineating spatially dependent objects, such as hydrocarbon accumulations from seismic data|
    US6882997B1|1999-08-25|2005-04-19|The Research Foundation Of Suny At Buffalo|Wavelet-based clustering method for managing spatial data in very large databases|
    DE19943325C2|1999-09-10|2001-12-13|Trappe Henning|Process for processing seismic measurement data with a neural network|
    US6295504B1|1999-10-25|2001-09-25|Halliburton Energy Services, Inc.|Multi-resolution graph-based clustering|
    US6226596B1|1999-10-27|2001-05-01|Marathon Oil Company|Method for analyzing and classifying three dimensional seismic information|
    US6574566B2|1999-12-27|2003-06-03|Conocophillips Company|Automated feature identification in data displays|
    FR2808336B1|2000-04-26|2002-06-07|Elf Exploration Prod|METHOD OF CHRONO-STRATIGRAPHIC INTERPRETATION OF A SEISMIC SECTION OR BLOCK|
    US6363327B1|2000-05-02|2002-03-26|Chroma Graphics, Inc.|Method and apparatus for extracting selected feature information and classifying heterogeneous regions of N-dimensional spatial data|
    US6618678B1|2000-05-26|2003-09-09|Jason Geosystems B.V.|Method of joint analysis and interpretation of the subsurface from multiple seismic derived layer property data sets|
    US6625541B1|2000-06-12|2003-09-23|Schlumberger Technology Corporation|Methods for downhole waveform tracking and sonic labeling|
    FR2813959B1|2000-09-11|2002-12-13|Inst Francais Du Petrole|METHOD FOR FACILITATING THE RECOGNITION OF OBJECTS, ESPECIALLY GEOLOGICAL, BY A TECHNIQUE OF DISCRIMINANT ANALYSIS|
    US6950786B1|2000-10-10|2005-09-27|Schlumberger Technology Corporation|Method and apparatus for generating a cross plot in attribute space from a plurality of attribute data sets and generating a class data set from the cross plot|
    US7006085B1|2000-10-30|2006-02-28|Magic Earth, Inc.|System and method for analyzing and imaging three-dimensional volume data sets|
    US6757614B2|2000-12-18|2004-06-29|Schlumberger Technology Corporation|Seismic signal processing method and apparatus for generating correlation spectral volumes to determine geologic features|
    US7203342B2|2001-03-07|2007-04-10|Schlumberger Technology Corporation|Image feature extraction|
    US6751558B2|2001-03-13|2004-06-15|Conoco Inc.|Method and process for prediction of subsurface fluid and rock pressures in the earth|
    US6473696B1|2001-03-13|2002-10-29|Conoco Inc.|Method and process for prediction of subsurface fluid and rock pressures in the earth|
    FR2824148B1|2001-04-30|2003-09-12|Inst Francais Du Petrole|METHOD FOR FACILITATING TRACKING OVER TIME OF THE DEVELOPMENT OF PHYSICAL STATES IN A SUBTERRANEAN FORMATION|
    DE10142785C2|2001-08-31|2003-07-03|Henning Trappe|Method for determining local similarity from 3D seismic measurement data|
    US6957146B1|2001-12-24|2005-10-18|Rdsp I, L.P.|System for utilizing seismic data to estimate subsurface lithology|
    FR2841344B1|2002-06-19|2005-04-29|Tsurf|METHOD, DEVICE AND PROGRAM PRODUCT FOR SMOOTHING SUBSURFACE PROPERTY|
    US7188092B2|2002-07-12|2007-03-06|Chroma Energy, Inc.|Pattern recognition template application applied to oil exploration and production|
    US7184991B1|2002-07-12|2007-02-27|Chroma Energy, Inc.|Pattern recognition applied to oil exploration and production|
    US20060184488A1|2002-07-12|2006-08-17|Chroma Energy, Inc.|Method and system for trace aligned and trace non-aligned pattern statistical calculation in seismic analysis|
    US7162463B1|2002-07-12|2007-01-09|Chroma Energy, Inc.|Pattern recognition template construction applied to oil exploration and production|
    US7308139B2|2002-07-12|2007-12-11|Chroma Energy, Inc.|Method, system, and apparatus for color representation of seismic data and associated measurements|
    US7295706B2|2002-07-12|2007-11-13|Chroma Group, Inc.|Pattern recognition applied to graphic imaging|
    US20050288863A1|2002-07-12|2005-12-29|Chroma Energy, Inc.|Method and system for utilizing string-length ratio in seismic analysis|
    GB2394050B|2002-10-07|2005-11-23| Westerngeco Seismic Holdings Limited|Processing seismic data|
    US7053131B2|2002-12-03|2006-05-30|Kimberly-Clark Worldwide, Inc.|Absorbent articles comprising supercritical fluid treated HIPE, I-HIPE foams and other foams|
    US7697373B1|2003-11-25|2010-04-13|Michael John Padgett|Method for deriving 3D output volumes using dip vector analysis|
    US7453767B1|2003-11-25|2008-11-18|Michael John Padgett|Method for deriving a 3D GRAZ seismic attribute file|
    US7453766B1|2003-11-25|2008-11-18|Michael John Padgett|Method for deriving 3D output volumes using summation along flat spot dip vectors|
    US7463552B1|2003-11-25|2008-12-09|Michael John Padgett|Method for deriving 3D output volumes using filters derived from flat spot direction vectors|
    US7206782B1|2003-01-29|2007-04-17|Michael John Padgett|Method for deriving a GrAZ seismic attribute file|
    US6754380B1|2003-02-14|2004-06-22|The University Of Chicago|Method of training massive training artificial neural networks for the detection of abnormalities in medical images|
    US7248539B2|2003-04-10|2007-07-24|Schlumberger Technology Corporation|Extrema classification|
    US6804609B1|2003-04-14|2004-10-12|Conocophillips Company|Property prediction using residual stepwise regression|
    US6970397B2|2003-07-09|2005-11-29|Gas Technology Institute|Determination of fluid properties of earth formations using stochastic inversion|
    GB0318827D0|2003-08-11|2003-09-10|Bg Intellectual Pty Ltd|Dip value in seismic images|
    US20080270033A1|2003-08-19|2008-10-30|Apex Spectral Technology, Inc.|Methods of hydrocarbon detection using spectral energy analysis|
    US7243029B2|2003-08-19|2007-07-10|Apex Spectral Technology, Inc.|Systems and methods of hydrocarbon detection using wavelet energy absorption analysis|
    US7092824B2|2003-10-20|2006-08-15|Ascend Geo Llp|Methods and systems for interactive investigation of geophysical data|
    US6941228B2|2003-12-05|2005-09-06|Schlumberger Technology Corporation|Method and system and program storage device for analyzing compressional 2D seismic data to identify zones of open natural fractures within rock formations|
    US20050171700A1|2004-01-30|2005-08-04|Chroma Energy, Inc.|Device and system for calculating 3D seismic classification features and process for geoprospecting material seams|
    US20060115145A1|2004-11-30|2006-06-01|Microsoft Corporation|Bayesian conditional random fields|
    WO2006100674A2|2005-03-21|2006-09-28|Yeda Research And Development Co. Ltd.|Detecting irregularities|
    FR2884636B1|2005-04-15|2007-07-06|Earth Resource Man Services Er|PR0CEDE OF HIERARCHICAL DETERMINATION OF COHERENT EVENTS IN AN IMAGE|
    US7266041B1|2005-06-21|2007-09-04|Michael John Padgett|Multi-attribute background relative scanning of 3D geophysical datasets for locally anomaluous data points|
    US20070067040A1|2005-09-02|2007-03-22|Anova Corporation|Methods and apparatus for reconstructing the anulus fibrosus|
    US7881501B2|2005-12-13|2011-02-01|Calgary Scientific Inc.|Local dominant wave-vector analysis of seismic data|
    EP2372398A3|2006-06-21|2012-04-11|Terraspark Geosciences, LLC|Extraction of depositional systems|
    CA2661449C|2006-08-31|2015-11-24|Shell Canada Limited|Interpreting a plurality of m-dimensional attribute vectors assigned to a plurality of locations in an n-dimensional interpretation space|
    DE102006050534B4|2006-10-26|2013-12-12|Airbus Operations Gmbh|Conduit system for an aircraft, in particular an aircraft|
    FR2909185B1|2006-11-27|2009-01-09|Inst Francais Du Petrole|METHOD OF STRATIGRAPHIC INTERPRETATION OF SEISMIC IMAGES|
    US9835743B2|2006-11-28|2017-12-05|Magnitude Spas|System and method for seismic pattern recognition|
    GB2444506C|2006-12-06|2010-01-06|Schlumberger Holdings|Processing of stratigraphic data|
    US8346695B2|2007-03-29|2013-01-01|Schlumberger Technology Corporation|System and method for multiple volume segmentation|
    US7658202B2|2007-05-08|2010-02-09|Kohler Co.|Low-profile valve assembly|
    EP2291686A4|2008-05-09|2017-01-04|Exxonmobil Upstream Research Company|Method for geophysical and stratigraphic interpretation using waveform anomalies|
    AU2009249473B2|2008-05-22|2014-01-16|Exxonmobil Upstream Research Company|Seismic horizon skeletonization|
    US8577613B2|2008-07-01|2013-11-05|Schlumberger Technology Corporation|Effective hydrocarbon reservoir exploration decision making|
    GB2465079B|2008-08-06|2011-01-12|Statoilhydro Asa|Geological modelling|
    AU2010257595B2|2009-06-09|2013-12-05|Shell Internationale Research Maatschappij B.V.|Method for analysing seismic data with attribute functions|
    US8463551B2|2009-11-17|2013-06-11|Schlumberger Technology Corporation|Consistent dip estimation for seismic imaging|
    US8326542B2|2009-11-19|2012-12-04|International Business Machines Corporation|Method and system for retrieving seismic data from a seismic section in bitmap format|
    US9607007B2|2009-12-23|2017-03-28|Schlumberger Technology Corporation|Processing of geological data|
    US8358561B2|2010-04-13|2013-01-22|Spectraseis Ag|Bayesian DHI for seismic data|
    US8380435B2|2010-05-06|2013-02-19|Exxonmobil Upstream Research Company|Windowed statistical analysis for anomaly detection in geophysical datasets|
    US9134443B2|2010-05-14|2015-09-15|Schlumberger Technology Corporation|Segment identification and classification using horizon structure|
    RU2573166C2|2010-05-28|2016-01-20|Эксонмобил Апстрим Рисерч Компани|Method for seismic analysis of hydrocarbon systems|
    US8515678B2|2010-07-29|2013-08-20|Schlumberger Technology Corporation|Chrono-stratigraphic and tectono-stratigraphic interpretation on seismic volumes|
    US8447525B2|2010-07-29|2013-05-21|Schlumberger Technology Corporation|Interactive structural restoration while interpreting seismic volumes for structure and stratigraphy|
    US20120090001A1|2010-10-07|2012-04-12|Tzu-Chiang Yen|Apparatus having multimedia interface and network access management integrated therein|
    GB201022128D0|2010-12-31|2011-02-02|Foster Findlay Ass Ltd|Active contour segmentation|
    US8861309B2|2011-01-31|2014-10-14|Chevron U.S.A. Inc.|Exploitation of self-consistency and differences between volume images and interpreted spatial/volumetric context|
    US8838391B2|2011-01-31|2014-09-16|Chevron U.S.A. Inc.|Extracting geologic information from multiple offset stacks and/or angle stacks|
    US20120197613A1|2011-01-31|2012-08-02|Chevron U.S.A. Inc.|Exploitation of self-consistency and differences between volume images and interpreted spatial/volumetric context|
    US8972195B2|2011-01-31|2015-03-03|Chevron U.S.A. Inc.|Extracting geologic information from multiple offset stacks and/or angle stacks|
    US9121968B2|2011-01-31|2015-09-01|Chevron U.S.A. Inc.|Extracting geologic information from multiple offset stacks and/or angle stacks|
    US20120322037A1|2011-06-19|2012-12-20|Adrienne Raglin|Anomaly Detection Educational Process|
    US8954303B2|2011-06-28|2015-02-10|Chevron U.S.A. Inc.|System and method for generating a geostatistical model of a geological volume of interest that is constrained by a process-based model of the geological volume of interest|
    US20130138350A1|2011-11-24|2013-05-30|Manoj Vallikkat Thachaparambil|Enhancedtermination identification function based on dip field generated from surface data|US10995592B2|2014-09-30|2021-05-04|Exxonmobil Upstream Research Company|Method and system for analyzing the uncertainty of subsurface model|
    US10267934B2|2015-01-13|2019-04-23|Chevron U.S.A. Inc.|System and method for generating a depositional sequence volume from seismic data|
    AU2017343749A1|2016-10-14|2019-04-11|Chevron U.S.A. Inc.|System and method for seismic facies identification using machine learning|
    WO2018080745A1|2016-10-28|2018-05-03|Exxonmobil Upstream Research Company|Method and system for enhancing meshes for a subsurface model|
    US11231515B2|2016-11-17|2022-01-25|Landmark Graphics Corporation|Multi-Z horizon interpretation and visualization from seismic data|
    CN109425889B|2017-08-22|2020-10-20|中国石油化工股份有限公司|Method for depicting ancient karst underground river|
    法律状态:
    2015-03-05| FGA| Letters patent sealed or granted (standard patent)|
    优先权:
    申请号 | 申请日 | 专利标题
    US201261598233P| true| 2012-02-13|2012-02-13||
    US61/598,233||2012-02-13||
    PCT/US2013/021179|WO2013122704A1|2012-02-13|2013-01-11|System and method for detection and classification of seismic terminations|
    [返回顶部]